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Abstract 

Using fermionic techniques, we compute exactly the large deviation function (ldf) 
of the time-integrated injected power in several one-dimensional dissipative systems of 
classical spins. The dynamics are T = Glauber dynamics supplemented by an injection 
mechanism, which is taken as a Poissonian flipping of one particular spin. We discuss the 
physical content of the results, specifically the influence of the rate of the Poisson process 
on the properties of the ldf. 

Keywords: Glauber dynamics, spin systems, large deviation functions, free fermions, 
stochastic systems, fluctuation theorem. 

1 Introduction 

Fluctuations in nonequilibrium systems have attracted interest of physicists and math- 
ematicians in recent years, due to fortuitous conjunctions of experiments [H [2], which 
aimed at measuring and understanding the fluctuations of global variables in hydrody- 
namic experiments, and theoretical works O [U [5], having stated new relations for the 
entropy production in some classes of nonequilibrium systems (the fluctuation theorems). 
This simultaneity was actually the cause of some bold assertions surmising a wider ap- 
plication of the so-called fluctuation theorems to dissipative systems [H [7] , whereas these 
relations were nothing but a manifestation of the time-reversal symmetry of the bulk 
dynamics. These conjectures were enforced by the apparent validity of the fluctuation 
relations in experiments and simulations, which yields incidentally a characteristic energy 
sometimes daringly termed "nonequilibrium temperature". Despite the fact that it has 
been convincingly proven that this apparent validity is only due to our inability of prob- 
ing the large deviation functions of injected power in zones of negative injected power 
[El El EH El [El H3], it is however interesting to show explicit examples of dissipative 
systems where the large deviation function (ldf) of injected power can be exactly com- 
puted, all the more so there exists very few such exact calculations on dissipative systems 
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[9j [TUl HH [151 [S] , for which the absence of detailed balance condition makes the situation 
more difficult to handle. 

Moreover, for conservative systems in nonequilibrium stationary states, interesting 
relations on the related ldf of integrated current have been discovered recently |17] . and 
it is interesting to compare the classes of conservative and dissipative systems as regards 
the fluctuations of the global variable associated with the external excitation (the very 
cause of the nonequilibrium state). 

In this paper, we consider a ID system of 2N (N — > oo) classical spins on a ring, la- 
belled from — N to N—l, updated according a T = Glauber dynamics, and supplemented 
by an injection mechanism. The T = Glauber dynamics is defined as follows : the prob- 
ability for a spin Sj to flip between t and t+dt is given by dt[l — Sj (sj_i + Sj+\)/2] , that is, 
Sj flips with a rate 1 if its neighbours are in different states, cannot flip if Sj = Sj-i = Sj+\ 
and flips with a rate 2 if — Sj = s 3 _x = s j+i- It is important to mention that these dynam- 
ics are dissipative: a flip either does not modify the energy of the system, or lessens it. 
As a result, energy must be injected by an additional mechanism of injection in order to 
drive the system into a non trivial stationary state. We will consider two different models 
of injection : for both models (labelled I and II hereafter), only the zeroth spin sq (the 
"external boundary" of the system) flips randomly as a Poisson process with a parameter 
A. The difference between the two models lies in the ability of the zeroth spin, besides the 
Poissonian flip, to spontaneously flip under the Glauber rules (model II) or not (model 
I); it is qualitatively more important than it could appear, as in one case (model II) both 
halves of the system are dynamically connected in the vicinity of so when in the other 
(model I) they are disconnected. In the thermodynamic limit, Model I can be viewed as 
two independent identical systems driven out of equilibrium by sq . 

For these systems, we computed exactly f(p), the large deviation function of the 
injected energy, a major observable associated with the fluctuations of the energy flux in 
stationary systems. It is defined as 

f(p) = lim r _1 log Prob((energy injected between and r)/r = p) (1) 

T—tOO 

As the dynamics are stochastic, no explicit formula can be put forward for II the energy 
injected in the system by A from t = to t = t, but fortunately such expressions are 
not needed for our purpose. In section [21 we expose the dynamics of the models in 
greater detail, as well as the computation and the exact results for g(a), the ldf of the 
characteristic function. In section [3] we numerically solve the inverse Legendre transforms 
and get the ldfs for the injected energy for the two models considered here. We discuss 
their physical properties, in particular their variations with respect to the change of the 
injection rate. 

2 Fermionic approach to the injected power 

It is useful to describe spin systems in the dual representation of the domain walls : be- 
tween the site and 1 is located the locus of a possible domain wall labelled 0, and so forth. 
The state of the system is thus characterized by C = (U-n, . . . , n_i, no, ni,ii2, • • • , Wjv-i), 
where the rii are either (no domain wall) or 1. There is 2 2N possible states in this rep- 
resentation, in contrast with the spin representation where a degeneracy s <-> —s doubles 
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this number. The dynamical equation for the probability is given by 



d t P(C) = X[P(C ) - P{C)\ + Y).P{Ci>{Cj - C) - P(C)w(C -» Cj)] (2) 

i 

where Cj holds for the state C whose domain wall variables rij and %_i have been changed 
(according to n — ► 1 — n). The T = Glauber dynamics corresponds to 

w(Cj — > C) = 2 — nj — nj_i (3) 
iu (C — > Cj) = nj + nj_i (4) 

where rtj and rij-i are the variables associated with the state C (we use this convention 
hereafter). Note that models I and II differ in the preceding equations in the way the 
summation over j is carried out: either it is not restricted (model II), or index j = is 
implicitely removed (model I). 

We consider that each domain wall contributes as a quantum of energy 1 to the global 
energy of the system. We are interested in the energy II injected into the system up to time 
t by the poissonian injection. Following [18] . the route to this time integrated observable 
begins with the consideration of the joint probability P(C,H,t), the probability for the 
system to be in the state C at time t having received the energy II from the injection. 
The dynamical equation for this quantity is readily 

d t P(C, IT) = \{P(C , n - 2)n n^ + P(C , n + 2)(1 - no)(l - n_i) 

+P(C ,n)[(l-n )n_i+(l-n-i)no]-P(C,n)}+J][P(C i ,n)«;(Cj ^ C)-P(C,n)w(C - C,-)] 

j 

(5) 

We define next the generating function of II as 

oo 

F(C)= e aU P(C,n) (6) 

n=-oo 

This quantity, summed up over the states, yields the generating function (exp(all)) from 
which one derives its ldf g(a): 

{e aU ) ~ e tg( - a) (7) 

This ldf g(a) is closely related to f(p), the ldf of the probability density function of IT, as 
they are Legendre transform of each other [9] : 

PTob(U/t = p) oc exp(*/(p)) (8) 

t— too 

f(p) = min (g(a) - up) (9) 

a 

Let us write the dynamical equation for F : 
d t F(C) = 

A [e 2a F(C )n n^ + e - 2Q F(C )(l - n )(l - n_i) + F(C )[(1 - n )n_! + (1 - n_x)n ] - F(C)] 

+ Y,[F(Cj)(2 - rij - - F(C)(nj + n^-i)] (10) 



3 



The function g(a) is closely related to the dynamical matrix at work in the rhs of the pre- 
ceding equation, since it is in general its largest eigenvalue. Thus, succeeding in computing 
g(a) is not as complicated as knowing the whole dynamics (eigenvalues and eigenvectors), 
but nevertheless it remains a difficult challenge, as the whole spectrum of the matrix has 
to be known. 

Fortunately, our problem belongs to the category of the "free-fermions" problems, for 
which a diagonalization of the dynamics into independent "modes" can be achieved. To 
do so, we construct the state vector as 

\<t>) = Y,nc)\c) ai) 

c 

where \C) is a vector in the 2 2N+l dimensional tensorial product of the space of states of 
domain walls (each of dimension 2). The dynamical evolution of \4>) yields expressions 
like F (Co)non-i\C) which can be represented as 

F(C )n n_ 1 |C) = F(Co)non-i\C) (12) 
c c 

= n n_ 1 ^F(C)|C ) (13) 
c 

= n n- 1 a^a 3 L 1 \<f>) = s^sZM (14) 

where 

«-(!!)■ '-(Si)- ---OS) <«> 

(we will have use also of s + = t (s~) ; note that s~s + = n). Note that for the time being, 
operators with different particle indices commute. Similar computations on other terms 
allows to write the dynamical equation for \<p) as 

d t \4>) = H\<f>) where (16) 

+ Y,i 2s U s t + s U s J + °7-i*t ~ s J s t ~ "J-i'U ^ 

3 

The fermionization of this "Hamiltonian" proceeds using the Jordan- Wigner transforma- 
tion 

C- N = a± N (18) 
def 

(19) 



where 



Cj = s+a z _ N a z _ N+1 . . . u]_ x (20) 

def J 

c] = sja z _ N a z _ N+1 . . . g z _ x (21) 



°' = \\ -I J < 22 > 
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is the third Pauli matrix. A major change provided by this transformation is that now 
operators with different indices do no longer commute. Instead, it can be verified that for 
all i and j 

{c t , Cj } = (23) 
{4,ct} = (24) 
{4,c j } = 6 itj (25) 

where the bracket is defined by {a, b} = ab+ba. These relations, characteristic of fermionic 
operators, yields a useful rewriting of the dynamical operator 

H = A[e 2o, cl 1 4 + e~ 2Q c c_i + cjc_i + c^co - 1] + ^[2cjCj_i + cjcj-i + c]-i c i ~ c ) c 3 ~ c )-i c i-^\ 

(26) 

The change of variables = e a c\c = e~ a c is compatible with the fermionic structure 
and allows the rewriting of the preceding equation (we omit the tildas in the following) : 

H = Afc^cJ + c c_i + cjc_i + ct^co - 1] + ^[2e 2Q CjCj_i + c]cj_i + c]_ 1 Cj - cjcj - ct_ 1 c J -_i] 

(27) 

In order to achieve the diagonalization, we follow the route described in the Appendix. 
The Hamiltonian can be written as 



e-X (28) 



where e is defined as: 

(model I) 



1 (model II) ( 29) 



The largest eigenvalue of H is thus given by 



{a) = -^\k q \-2N+l-£-\ (30) 



9 - 2 



provided all the A q are real. If one takes into account possible complex eigenvalues, g(a) 
will be rather defined by 



5(«) = ^|Re(A g )|-2iV + l- e -A (31) 
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As explained in the appendix, {A q , —A q } q is the spectrum of the matrix 

M »=(o -a) (32) 

with the notations of the appendix; moreover we can compute its characteristic polynomial 
from ([93]) . taking advantage of B almost empty. We define E{n) = (A + ^l)~ 1 D{A — 
and (f93l) shows that we have to compute det(l + BE); the matrix (1 + BE) is quite 
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simple, since B = A(M_i o — Mq _i), where {Mij}ij is the canonical basis of the 2N x 2N 
matrices. As a result BE is zero except on the lines —1 and 0, and one has 



det(l + BE) = (1 - XE_i, )(l + \E -i) + A 2 £ ,o#-i,- 



(33) 



= (1 + A^o,-i(-/^))(l + A£ ,-iGu)) + A 2 £ ,o(^)#-i,-i(aO (34) 

where we exploited the fact that E T (/j,) = —E(—n). Note in passing that the symmetry 
of this expression with respect to [i — > —fi, demonstrated with general arguments in the 
appendix, is blatant here, as the Ejj are antisymmetric functions of \i. 

A further simplification is provided by the fact that the physical system has itself a 
symmetry i *-* — i — 1 (expressed in the matrices A, B and D as Xij = which 
gives E 00 (n) = -E-i -i(fi) and £7 ,-i(m) = #0.-1 (-/■*), whence 



det(l + BE) = [1 + A£ ,-iUO + A£7 ,o(/i)][l + Af^-iM - XEqM] 



(35) 



The next step is to compute the elements (0, j) and (— 1, j) of the matrix (/ild — A)^ 1 , 
or equivalently the associated cofactors of fild—A, termed Cof(i, j). Precisely, we need the 
cofactors Cof(0, j), and the others are obtained using Cof(i, j) = Cof(j, i) and Cof(i, j) = 
Cof(-t- 1). We get : 



Cof(0,j) ^ 
det(/i - A)n^ 



[\ + e]{-l)-ixj +1 

(x + e- l) 2 - (A + e) 2 
[x + e- l](-lpx-^ 

I (x + e-l) 2 -(A + e) 2 



if J < 
if J ^0 



(36) 



where 



x 



+ 2 ± yV + 4/i 



(± such that |x| be maximum) 



After cumbersome computations, we get, for N ^ oo : 

(x+x_) 2iV 



Xw (m) = G+(fi)G-(fi)- 4 

A 1 



16^ 2 



= (x+ - 1 + 2s + A)(x_ - 1 - A) + A 
GL(/z) = (x_ - 1 + 2e + A)(x+ - 1 - A) + A 



A + 2ee 2a + 2e 



'2a 



X_ 



x + 



A + 2ee 2a + 2e 



x+x_ — 1 

2a x + ~ x ~ 



X_i_X_ 



1 



X+ 



X- 



/i + 2 + V/x 2 + 4// 



for Re(/i) > -2 



-»+ 2 -f^ forRe(/i)>2 



(37) 

(38) 
(39) 
(40) 

(41) 

(42) 



Finally, one has to combine this last result with (|3ip . This is done via the formula 
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2iV + 1 — e — A 



(43) 



where the contour for the complex integration is a big (big enough to enclose all singular- 
ities of the meromorphic function) half-circle leant on the imaginary axis, with its belly in 
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the Re(//) > region, followed counterclockwise. The logarithmic derivative of xm yields 
different additive contributions from formula (|38|) . Let us consider first the contribution 
(x+X-) 2N . The integral 

2N 



4iir 



j> dfj, 


x' + x' 


y, — 







(44) 



can be exactly computed, as the analytical singularities in the right half plane of x+ and 
x_ are located only on the real segments [—4, 0] and [0, 4] respectively; as a result, we can 
compute this integral with a contour encircling [0, 4] and sticking to it. Yet some care 
must be taken however, for the actual expression of X- is not uniform on the whole half 
plane, and notably on the real axis: whereas x+ is everywhere given on Re(/i) > by 
{fi + 2 + y V + 4/i]/ 2, x_ is given by [-/j, + 2 + y^ 2 - 4/j]/2 for Re(fi) G [0, 2] and by 
[— fi + 2 — \J ii 2 — 4:fj]/2 for Ke(fi) > 2. With all these precautions we get the simple result 
I = 2N. This extensive term vanishes with the corresponding term coming from Tr(A), 
which is coherent with the injection properties being asymptotically independent of the 
size of the system (in our ID model) . 

We can also compute easily the contribution of the term (/x 4 — lftfi 2 )^ 1 . The associated 
integral yields —2, thence we can recast the expression for g(a) into 



G+(/x) G-Qj) 



1 — e — A (45) 



2.1 The case a = 



It is possible to verify partially the correctness of the formula. Indeed, when a = 0, we 
know the value of g(a): g(0) = 0, that maximum eigenvalue being associated with the 
stationary distribution (encoded in the corresponding eigenvector). Thus, we must verify 
this expected result, which holds whatever the values of A and e. 

To do that, we rewrite the expressions for G+ and G- in a more transparent form. 
Using the relations x\ — (±/x + 2)x± + 1 = 0, one can write 

e 2a - 1 

G+(n) = (x- - 1) [(x+ - 1)(1 + 2A/V) + 2e) + 2A T [e(x+x. - l) + x _ - x+] 

X+X- — 1 

(46) 

G- is easily obtained by switching x + and x_, and making [i — * — 

e 2a - 1 

G_(/x) = (x+ - 1) [(z_ - 1)(1 - 2A///) + 2e] + 2A -[e^+x. - 1) +x+ - x_] 

x+x_ — 1 



(47) 



2.1.1 £ = 0,ct = 



This case is the simplest, since G± = (x_ — l)(x+ — 1)(1 ± 2A///) and the logarithmic 
derivative yields six contributions. The two contributions coming from {x + — 1) and the 
contribution coming from (1 + 2A//i) yield zero, as no branch cut/poles are associated 
with these functions in the half plane ReQii) > 0. Thus, from (|45p . 



^^-^ + -L£)-i-a 



fm —} fll,ll {x TV) 1 A 

1 /" 2x' 

k-J^—- 1 < 49 » 
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The contour of integration can be stuck to the branch cut [0,4] (that is [0^4] — iO 
followed by [4 -> Q]+ iQ. From - iO to 2 - iO, x_(/z = y-i0) = [-^ + 2 + ^2-4^/2 = 
\—y + 2 + i\J Ay — y 2 ]/2; from 2 — iO to 4 — iO, the expression for x_(/z) is changed, but 
the limit for fx sticking below the branch cut is such that one has still x-(fi = y — iO) = 
\—y + 2 + i\J Ay — y 2 ]/2. The expressions above the branch cut are just the complex 
conjugates. As a result x_ covers counterclockwise the set of complexes of modulus 1, 
from z = 1 (fi = — iO) to z = — 1 (/z = 4 — iO) with positive imaginary parts, and from 
2 = — 1 (n = 4 + iO) back to z = 1 with negative imaginary parts. As a result, we get the 
expected value: 



1 



5(0) = ^~ 



2m 



dz- 



z — z 



1*1=1 



z- 1 







(50) 



2.1.2 e = l,a = 

The computation in that case is slightly different, in particular as regards the localization 
of the singularities. The terms (x+ — 1) and (x+ + l + 2A^~ 1 (x + — 1)) have no singularities 
in the half plane under consideration (the singularity \i = is avoided). The term (x_ — 1) 
yields a factor 1/2, as already computed. The term 



x_ + 1- 2A^ x (x_ - 1) 



x_ — 1 



(x 2 + 2Ax_ - 1) 



(51) 



is the richest, as it has not only a branch cut, but also displays an isolated pole for 
H = 2[l + VA 2 + 1]. Thus, 



5(0) = -4^f^ + if^x: 



dz (z - l) 2 (z + A) 
1 2ivr /| Z | =1 T z 2 + 2Az - 1 



- 1 + Va 2 + 1 - A 
(52) 

(53) 



= -J-/ ^ (izi^ - 1 + 0^1- A = O 
2z7r 7| z |=i z z + 2Az — 1 

2.2 The general case (a 7^ 0) 

The preceding particular cases show clearly that the analysis of the singularities becomes 
very complicated in the general the new term proportional to e 2a — 1 mixes x+ 

and X-. It is thus more appropriate to switch to an integral representation of the result. 
To achieve this in a simple way, we make use of the preceding verifications. We define 



2A(e 2Q -l) l + e(x-x+-l)/(x--x+) 



x-x + + 1 fi/2 + A + £fj,/(x+ — 1) 

and similarly L_ using the rules x_ <-> x + and <-> —/it. Using the identity 

(x_ — x + )(x_x + + 1) 



(x+ — l)(x_ — l)(x_x + — 1) 



2/fi 



one shows that 



1 

4i7T 



+ 

U I- 



(54) 



(55) 



(56) 
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2.2.1 Model I (e = 0) 

The simplest expression is provided by the system I (e = 0) since in that case: 

i 2A(e 2a -l) 1 
I± = 1 + — x 



x_x+ + l ±n/2 + X 



(57) 



The integrand of (|56|) is of order /U -3 , so the half circle of integration can be made infinite, 
the only non vanishing contribution being that of the vertical axis Re(/i) = 0. It is easy 
to verify that on this axis x + = x*_ and X-(-iy) = x_(iy)*. As a result, one has 



9(a) 



1 

Air 



dy y 



DO 

oo 



1 dU 1 dl- 

+ 



1+ dy I- dy 



The explicit expression for |/+| 2 is: 



|/ + (iy)| 2 = l + 



and 



2a\ 



AX 2 (e 2a - l)(x_x+ + e 
(x_x+ + 1) 2 (A 2 + y 2 /4) 



x_x+ = |x + | 2 = - + 2 + v 7 -?/ 2 + Aiy 

2 + ^-y + ^y 2 + 16 ) + (y + -^|^ + ^ 2 + 16 ) 

To summarize, we can write the result as 



1 

y^O 4 



2 /°° , . / A 2 (e 2 "-l)Wu) + e 2 ») 
7r7 V + l) 2 (A 2 /4 + n i ) 



1 + V2\J-u 2 + yV + m 2 ^) + (lu + u 2 + Vu 4 + u 2 ^j 



A useful relation for ip(u) can be deduced: from 
1 1 + x + /x_ 



X_X + + 1 (1 + x_x + )(l + x + /x- 



1 1 
— + — 

X + X- 



one gets 



1 1 



ij>(u) + 1 2 ^2 



i=\J -u 2 + \/w 4 + 



(58) 
(59) 



(60) 

(61) 
(62) 



(63) 
(64) 

(65) 
(66) 



This suggests the change of variables u = 1/sinhr. This change gives the alternate 
expression: 



2 coshr, /_ A 2 (e 2a -l)(coth 2 f + e 2 ") 



g(a) = - I dr 2 
vr ,/ smh r 



log 1 + 



(coth 2 ^ + l) 2 (^sinh 2 r + l 



sinh r 



(67) 
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A partial verification of the result can be done, via the computation of g'(0). This quantity 
equals the mean injected power in the stationary state, which can be computed directly 
by other means [20j. We have readily from (|63p : 

4^2 roc ^ u 

m = —Jo W«) + 1)(AV4 + ^) (68) 
= 2A(1 + A - VA 2 + 2A) (69) 

which is the expected result. Note that it is not surprising to find exactly the same result 
as in [20] despite the fact that here the system is duplicated with respect to the one studied 
in [2D]. There, the energy associated with a domain wall was twice the value adopted here 
(1 per domain wall). 



2.2.2 Model II (e = 1) 



This model yields an expression slightly more complicated for g(a), due to an involved 
expression for I + (formula (|54p ). Some algebraic manipulations allow only a partial sim- 
plification : 



2A(e 



2a 



where ip = This yields the result 



1) 2/# - (V - 1) 2 

X 
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du log 



1 + 



A(e 



2a 



vV 2 + 4^ + 2A 



1) l + i(if>- l) 2 /8uip 

X 



(70) 



^ + 1 X/2 + V-u 2 + iu 



(71) 



with -0 = V( u ) given by (1641) . This seemingly irreducible complexity we could have 
expected, since in contrast with model I, the mean injected power is here not computable 
by elementary operations (i.e. the dynamical equations for the correlators (crnCi) are not 
closed) . 



3 The large deviation functions 

In this section, we discuss the properties of the ldf f(p) associated with the integrated 
injected energy : 

Prob(n/i = p) oc exp(t/(p)) (72) 
f(p) = min (g(a) - ap) (73) 

a 

The typical shapes of the large deviation functions are plotted on figure [TJ Their 
abscissae are rescaled to align their maximum on 1. The missing information on (p) is 
plotted on figure [2] for both models (obtained via (p) = g'(a = 0)). The mean injected 
power behaves monotonously in both models (fig. [2]), with but a maximal asymptotic 
value twice larger (2) for model II. We verified the consistency of this result by a direct 
simulation of the model II (results not shown; for model I, we have the exact formula 
(p) = 2A(A + 1 — \/A 2 + 2A)). That model II induces an injection larger than model I 
is obvious: the spin sq can switch back spontaneously and swiftly just after the creation 
of two domain walls, putting the system back in a state ready to accept anew a positive 
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(a) model I (b) model II 




p/<p> p/<P> 

Figure 1: Typical shapes of the large deviation functions of model I (left) and II (right), for 
various values of A. The abscissae only are rescaled. For model II, A = oo is not plotted due to 
the particular analytical behaviour of / in that case (see text); however we believe A = 10 to 
be close to the asymptotic case. 




X 



Figure 2: Mean injected power as a function of A for the two models considered. 

injection of energy; in model I, the domain walls have to move away diffusively before 
such a state can be reached. It is interesting to note that for model II lim^oo g'(0) is not 
given by the simple limit taken inside the integral (which would have given (p) = 1). We 
postpone the discussion on this peculiarity to the end of the section. 

The ldfs for both models have the interesting property that they have a finite limit 
when p — > 0: arbitrary small integrated injected energies "are not so rare"; this feature 
would a priori favour the presence of a negative tail in the case where one averages over 
different initial conditions (as described in [9]); whether it is or not the case here is an 
issue beyond the scope of this paper. 

The ldfs are apparently close to a parabola. This is of course not true, and the 
deviation from the parabola is a rather important feature to look at. To this end, one 
defines <r(A) = \ f"({p})\ x (p) 2 , which inverse is a measure of the relative fluctuations of II 
(up to a factor). We see from figure[3]that <r(A) increases gently with A and reaches 

a constant value. This is consistent with the shrinking of / with increasing A to be seen 
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X 

Figure 3: Parameter cr(A) = \f"((p))\ X (p) 2 as a function of A. 
on figure [TJ 

We then use c(A) to rescale the ldfs so that the curvatures are equal to 1, namely 
consider f{p)/a{\) as a function of p/(p)- The results are plotted for model I on figure 
HI together with the reference parabola — \{p/{p) — l) 2 . Corresponding curves for model 
II are very similar, as can be seen on figure [5l where the ldfs for both systems with 
A = 2 (maximum discrepancy of parameters x, see later) are simultaneously plotted. 




p/<p> 

Figure 4: Rescaled large deviation functions for the model I. Note that the curves for A = 0.1 
and A = oo are almost identical, and this also true for A = 1 and A = 5. This is related to the 
non-monotonous behaviour of the parameter x(^)> see text- 

We remark that the ldf are rather different from a parabola, and display a marked tilt 
counterclockwise. Interestingly enough, the magnitude of the tilt is not a monotonous 
function of A. This can be simply caught by inspecting x(A) = g"'(0)g'(0)/g" 2 (0), the 
third coefficient of the Taylor expansion of the scaled ldf near the maximum: f{p)/a(X) + 
\{p/{p) — l) 2 — \x{p/{p) ~ I) 3 - As shown on figure x(A) has a non-monotonous 
behaviour with respect to A, and mainly dictates the location of the asymptotic tail. 
Note that this parameter is positive, in accordance with the counterclockwise tilt already 
mentioned. The physical origin of these features is rather complicated. The tilt is mainly 
associated to the temporal correlations of the process, which are revealed in the statistics 
of rare events only. The probability of an extra injection of energy is more likely than the 
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Figure 5: Rescaled large deviation functions for A = 2 in both models. 
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Figure 6: Asymmetry parameter x{ty- See text for details 

opposite, because of a kind of positive feedback : if for a while, an extra amount of energy 
has been already transfered into the system, the typical density of domain walls near 
the boundary experiences a positive fluctuation; as a result, subsequent incoming domain 
walls are more likely to eventually annihilate with an alter ego, instead of returning to the 
boundary (conversely a negative fluctuation is hard to maintain due to a growing shortage 
of domain walls to react with, see [9] for similar arguments). 

That this positive feedback has a minimum is related to the structure of the stationary 
state. On figure[7]is plotted the repartition function of the domain wall nearest to sq in the 
stationary state for model I (we restrict the discussion to model I for sake of simplicity) and 
various values of A — these pdfs were obtained numerically although the exact computation 
is a priori tractable. It is seen that the probability for the first domain wall to be between 
spin 1 and spin 2 (that is, next to the first possible domain wall, which so to speak does 
not belong to the bulk as it is directly affected by the boundary flips) is maximum for 
A ~ 1 — 2 (and is maximum with respect to the subsequent locations). This property 
is also true if one disregards the domain walls located in contact with the spin so and 
restricts oneself to "inner" domain walls (see inset of fig [7]). The correspondence between 
x(A) and the rescaled height of the peak 1-2 can be made quantitative by a correlation 
plot (fig. [3(b)), where a linear relation between the two quantities can be approximately 
drawn. 

Thus, a density fluctuation is typically localized very near the boundary especially 
for A ~ 2, and this proximity works against the correlations and the positive feedback 





2 4 6 8 10 
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(a) (b) 




spin index Peak 1 -2 rescaled 

Figure 7: (a) Repartition function of the domain wall nearest to the spin for model I. The 
inset shows the distribution function restricted (and rescaled in accordance) to the bulk domain 
walls, (b) Correlation plot between the height of the first peak (position 1-2) of the rescaled 
repartition function and xW- The line is a linear fit. 



effect explained above, as after just one move these domain walls can disappear through 
a negative energy injection event; this explains roughly the minimum of x(^) near A ~ 2. 
Finally one explains also the less pronounced effect for model II by the fact that the 
two halves of the system see each other in that case and domain walls from one side can 
pervade the other: this permeability probably reduces the correlations between density 
fluctuations and injection, thereby decreasing the asymmetry parameter. 




Figure 8: Nonequilibrium temperature T curv (A) = g"(0)/2g'(0). See text for details. 

Interestingly enough, this parameter \ has a similar behaviour as the nonequilibrium 
"temperature" one can define for such dissipative NESS from f(p) [TUl E] : T curv = 
g" (0)/2</ (0). Figure [8] shows this quantity for both models, and it is clear that this char- 
acteristic energy constructed from the fluctuations of II is sensitive to both the temporal 
correlations of the process (as x) an d the averages quantities like the mean injected power. 
As such, it bears a rather complicated physical content, in contrast with an ordinary tem- 
perature concept. 

Before to conclude, it is worth mentioning a mathematical subtlety of model II: the 
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naive A — ► oo limit yields 



2 f°° 
g(a) — — I du log 
Jo 



_ 2(e 2a -l)(l + i(^-l) 2 /8n^) 



V' + i 



whence one would have 



t lio V> + 1 

= 1 



(74) 

(75) 
(76) 



But this is not the correct result ({p) = 2), and it is easy to see why : the exact formula 
for finite A is, provided \/—u 2 + iu = R(u) + il(u) : 



8A 



00 , A/2 + R + IN) - l) 2 /8™/> 

fin 

,T ./o (^ + i)([A/2 + J R]2+/2) 



y^-u 2 + vV + 



I(u) = -^V^ 2 + V^T 
v2 



(77) 
(78) 
(79) 



The missing part comes from the third term of the numerator : / ~ u for large u, hence 
the associated integral is of order 1/A and not 1/A 2 as naively expected. This problem is 
ultimately due to the fact that an isolated pole exists for model II at the right hand side 
of equation (|56p (see section [2]), an analytical property extremely sensitive to the form 
of the approximants: some sub-dominants terms are of primarily importance to keep the 
analytical structure untouched and makes allowance for the asymptotic locations of the 
singularities. 



4 Conclusion 

In this paper, we were able to compute exactly by fermionic techniques large deviation 
functions of the injected power for two ID models of classical spins driven by a T = 
Glauber dynamics in the bulk and a Poissonian flip on the boundary. The results 
highlight the influence of the flipping rate on the low-frequency fluctuation properties of 
the injection; we showed in particular that the third derivative of / is a sensible measure of 
the interplay between on one hand density fluctuations and on the other hand the ability 
for the external operator to transfer efficiently energy inside the system. By the way, 
we found and explained why this third derivative (or equivalently the third cumulant of 
the time-integrated injected energy) is always positive, a property surprisingly found also 
in conservative diffusive systems for the integrated current (see formula (4) in [17]): this 
similar behaviour is a bit surprising, as inner dynamics of dissipative and conservative 
systems are by nature much different. 
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5 Appendix 



We give some details on the diagonalization of a Hamiltonian of the type: 



(80) 



where A is symmetric real and B and D antisymmetric real. We postulate the transfor- 
mation 



,t 



(81) 

c = U*$ + V*C (82) 

where U and V are 2N x 2N matrices. The fermionic structure is preserved if 

1 (83) 

(84) 



UU* T + VV* T 
UV T + VU T = 



It must be stressed that these equations are only 2N{2N + \) independent relations among 
2(2iV) 2 coefficients. Thus 27V(2iV - 1) relations can still be imposed to the coefficients of 
U and V. 

Next, we can replace (|8ip into (I80p and extract four different contributions : a contri- 
bution made with terms (N(2N — 1) terms) and £'£1 (_/V(2iV — 1) terms) we want to 
get rid of : 



-(£ f ) T [V T AU* - U* T AV + V T BV + U* T DU*} (tf) 
~(C) T [U T AV* - V* T AU + U T BU + V* T DV*] (£) 



(85) 
(86) 



this gives precisely the 2N{2N — 1) remaining relations on U and V (the matrices are 
antisymmetric); thirdly a contribution £T£ 



« t ) T M(0 = (£ t ) T 

del 



V T AF* - £/ T AET + ^{V T BU - U T BV) + ^(U* T DV* - V* T DU* 



(0 



(87) 



Our goal is to diagonalize this matrix. If this is possible, the diagonalization is completed 
using M = Q~ 1 AmQ and £ = Q£,C^ = Q~ But this is apparently not always 
possible, as the inspection of M clearly shows. Two special situations allow however a 
diagonalization: if B = —D, M is hermitian and therefore diagonalizable ; if U and V 
can be found real, M is symmetric and also diagonalizable. But the last case we cannot 
immediately recognize. In the following we assume that a diagonalization of M can be 
completed. 

There is a fourth term in the decomposition, which is a constant arisen from the 
commutation — > This constant reads 



Tr 



u t au* + -c/ T 5y + -v* t dit 



~Tr(M) + ^Tr(A) 



where relations among [7 and F were used. As a result, the Hamiltonian can be cast in 
the following form : 



(89) 
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where the A q are the eigenvalues of M. This formula differs slightly from those presented 
in [19] where the Tr(^4) term is not written. It is interesting to see that the precise sign 
of A q is irrelevant, as far as the spectrum of H is concerned. The eigenvalues of H are 
actually given by 



(90) 



where the e q are ±1. 



The preceding development essentially yields the conclusion embodied by the equation 
(|89p . provided that the diagonalization is possible. To study the diagonalization itself, it 
is particularly clever to follow the route described in [19J, where in (|81j) . the roles of c and 
£ are reversed. Starting from (f89j) and using [H, £ 9 ] = A q ^ q , one gets the following result: 
the spectrum of the matrix 

m "={d -a) < 91 > 

is made with the eigenvalues of M. Actually it is exactly the spectrum of M plus its 
symmetric with respect to zero : multiplying fjld — Mq by 



(fi + A^Difi-A)- 1 (ji + A) 



-in... „W ,,. , ^-l I (92) 



it is easily seen that the characteristic polynomial XAfbC/f) = det(/ild — Mo) can be recast 
into 

XmM = XA(riXA(-n)(-l) Aim{A) det(l - B(» + A)- l D{^ - A)" 1 ) (93) 

Next, using the relation det(l+AB) = det(l+BA) valid whatever A and B, and det(A) = 
det(A T ), it is readily seen that XM {fJ-) = Xm (— A 4 )- Thus, the spectrum of H is just 
duplicated in Mq according to Sp(Mo) = (Sp(-fT), — Sp(if)} 

The question of the diagonalization of H is closely related to the diagonalization of 
Mq. The case D = —B is explicit in this representation since in that case Mo is symmetric 
and thus diagonalizable with real eigenvalues. Actually, as the non diagonalization of a 
matrix is an accident rather than the rule, we can consider these cases as exceptional. 
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